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ABSTRACT 

We consider the variational problem of maximizing the weighted equilibrium 
Green's energy of a distribution of charges free to move in a subset of the upper 
half-plane, under a particular external field. We show that this problem admits a 
solution and that, under some conditions, this solution is an S-curve (in the sense of 
Gonchar- Rakhmanov) . The above problem appears in the theory of the semiclas- 
sical limit of the integrable focusing nonlinear Schrodinger equation. In particular, 
its solution provides a justification of a crucial step in the asymptotic theory of 
nonlinear steepest descent for the inverse scattering problem of the associated lin- 
ear non-self- adjoint Zakharov-Shabat operator and the equivalent Riemann-Hilbert 
factorization problem. 
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1. INTRODUCTION 

Let H = {z : Imz > 0} be the complex upper-half plane and H = {z : Imz > 
0} U {oo} be the closure of H. Let also K = {z : Imz > 0} \ {z : Rez = 0, < 
Imz < A}, where A is a positive constant. In the closure of this space, K, we 
consider the points ix + and ix^., where < x < A as distinct. In other words, we 
cut a slit in the upper half-plane along the segment (0, iA) and distinguish between 
the two sides of the slit. The point infinity belongs to K, but not K. We define F 
to be the set of all "continua" F in K (i.e. connected compact sets) containing the 
distinguished points + ,0_. 

Next, let p°(z) be a given complex-valued function on H satisfying 

p°(z) is holomorphic in H, 
p (z) is continuous in H, 

(1) 

Re[p°(z)] = 0, for z e [0,iA], 
Im[p°(z)] > 0, for z e [0,iA] Ul. 
Define G(z; n) to be the Green's function for the upper half-plane 

(2) G[z . v) ^i og \l^l 

\z - n\ 

and let du°(n) be the nonnegative measure —p°(n)dn on the segment [0, iA] oriented 
from to iA. The star denotes complex conjugation. Let the "external field" <j) be 
defined by 

(3) <f){z) = - J G(z;n)dpP{ri) - Re(inJ p°{n)dn + 2i,I(zx + zH)), 

where x,t are real parameters with t > and J = 1, for x > 0, while J = 
— 1, for x < 0. Re denotes the real part. 

The particular form of this field is dictated by the particular application to the 
dynamical system we are interested in. The conditions (1) are natural in view of 
this application. But many of our results in this paper are valid if the term zx + z 2 t 
is replaced by any polynomial in z. Here x, t are in fact the space and time variables 
for the associated PDE problem (see (9)-(10) below). 
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Let M be the set of all positive Borel measures on K, such that both the free 
energy 

(4) E(p)= J J G(x,y)dn(x)dfi(y), ^ei 
and / 4>d/j are finite. Also, let 

(5) V"(z) = J G{z,x)d f i(x), (iel 

be the Green's potential of the measure /i. 
The weighted energy of the field <p is 

(6) E+in) = E(ji) + 2 J 4>dn, (ieM. 

Now, given any continuum F G F, the equilibrium measure X F supported in F 
is defined by 

(7) E rj> (X F ) = mm^ eM(F )^(/x), 

where M(F) is the set of measures in M which are supported in F, provided such 
a measure exists. E<p(\ F ) is the equilibrium energy of F. 

The aim of this paper is to prove the existence of a so-called S-curve ([1]) joining 
the points 0+ and 0_ and lying entirely in K, at least under some extra assumptions. 
By S-curve we mean an oriented curve F such that the equilibrium measure \ F 
exists, its support consists of a finite union of analytic arcs and at any interior 
point of supp^i 

where the two derivatives above denote the normal (to suppfi) derivatives. 

To prove the existence of the S-curve we will first need to prove the existence of 
a continuum F maximizing the equilibrium energy over F. Then we will show that 
the maximizer is in fact an S-curve. 

It is not always true that an equilibrium measure exists for a given continuum. 
The Gauss-Frostman theorem ([2], p. 135) guarantees the existence of the equilib- 
rium measure when F does not touch the boundary of the domain H. This is not 
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the case here. Still, as we show in the next section, in the particular case of our 
special external field, for any given x,t and for a large class of continua F not 
containing infinity, the weighted energy is bounded below and X F exists. So, in 
particular, we do know that the supremum of the equilibrium weighted energies 
over all continua is greater than — oo. 

S-curves were first defined in [1] , where the concept first arose in connection with 
the problem of rational approximation of analytic functions. Our own motivation 
comes form a seemingly completely different problem, which is the analysis of the 
so-called semiclassical asymptotics for the focusing nonlinear Schrodingcr equation. 
More precisely, we are interested in studying the behavior of solutions of 

thd t ip + %d 2 ^ + M 2 V = o, 

(9) L 

under ip(x, 0) = ipo(x), 

in the so-called semiclassical limit, i.e. as H —> 0. For a concrete discussion, let us 
here assume that ipo (x) is a positive " bell-shaped" function; in other words assume 
that 

ip (x) > 0, iel, 
i>o(-x) = "4>o(x), 

(10) ipo has one single local maximum at 0, V'o(O) = A, 

K(0) < 0, 
ipo is Schwartz. 

This is a completely integrable partial differential equation and can be solved 
via the method of inverse scattering. The semiclassical limit is analyzed in the 
recent research monograph [3] . In Chapter 8 of [3] it is noted that the semiclassical 
problem is related and can be reduced to a particular "electrostatic" variational 
problem of maximizing the equilibrium energy of a distribution of charges that are 
free to move under a given external electrostatic field (assuming that the WKB- 
approximated density of the eigenvalues admits a holomorphic extension in the 
upper half-plane). In fact, it is pointed out that the existence and regularity of an 
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S-curve implies the existence of the so-called "g- function" necessary to justify the 
otherwise rigorous methods employed in [3]. 

We would like to point out that the problem of the existence of the " g- function" 
for the semiclassical nonlinear Schrodinger problem is not a mere technicality of 
isolated interest. Rather, it is an instance of a crucial element in the asymptotic 
theory of Riemann-Hilbcrt problem factorizations associated to integrable systems. 
This asymptotic method has been made rigorous and systematic in [4] where in 
fact the term " nonlinear steepest descent method" was first employed to stress the 
relation with the classical "steepest descent method" initiated by Riemann in the 
study of exponential integrals with a large phase parameter. Such exponential in- 
tegrals appear in the solution of Cauchy problems for linear evolution equations, 
when one employs the method of Fourier transforms. In the case of nonlinear 
integrable equations, on the other hand, the nonlinear analog of the Fourier trans- 
form is the scattering transform and the inverse problem is now a Riemann-Hilbert 
factorization problem. While in the "linear steepest descent method" the con- 
tour of integration must be deformed to a union of contours of " steepest descent" 
which will make the explicit integration of the integral possible, in the case of 
the " nonlinear steepest descent method" one deforms the original Riemann-Hilbert 
factorization contour to appropriate steepest descent contours where the resulting 
Riemann-Hilbert problems are explicitly solvable. 

In the linear case, if the phase and the critical points of the phase are real it 
may not be necessary to deform the integration contour. One has rather a Laplace 
integral problem on the contour given. For Riemann-Hilbert problems the analog 
is the self-adjointness of the underlying Lax operator. In this case the spectrum of 
the associated linear Lax operator is real and the original Riemann-Hilbcrt contour 
is real. The "deformation contour" must then stay near the real line. One novelty 
of the semiclassical problem for (9)-(10) studied in [3] however is that, due to the 
non-self-adjointness of the underlying Lax operator, the "target contour" is very 
specific (if not unique) and by no means obvious. It is best characterized via the 
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solution of a maximin energy problem, in fact it is an S-curve. The term "nonlinear 
steepest descent method" thus acquires full meaning in the non-self-adjoint case. 

Given the importance and the recent popularity of the " steepest descent method" 
and the various different applications to such topics as soliton theory, orthogonal 
polynomials, solvable models in statistical mechanics, random matrices, combina- 
torics and representation theory, we believe that the present work offers an im- 
portant contribution. In particular we expect that the results of this paper may 
be useful in the treatment of Riemann-Hilbcrt problems arising in the analysis of 
general complex or normal random matrices. 

On the other hand, we believe that the main results of this paper, Theorems 
3, 4, 5, 7, 8 are interesting on their own. This paper can be read without the 
applications to dynamical systems in mind. It concerns existence and regularity of 
a solution to an energy variational maximin problem in the complex plane. 

The method used to prove the existence of the S-curves arising in the solution of 
the " max-min" energy problem was first outlined in [1] and further developed in [5] , 
at least for logarithmic potentials. But, the concrete particular problem addressed 
in this paper involves additional technical issues. 

The main points of the proof of our results are: 

(i) Appropriate definition of the underlying space of continua (connected compact 
sets) and its topology. This ensures the compactness of our space of continua which 
is crucial in proving the existence of an energy maximizing clement. 

(ii) Proof of the semicontinuity of the energy functional that takes a continuum 
to the energy of its associated equilibrium measure (Theorem 3). 

(iii) Proof of existence of an energy maximizing continuum (Theorem 4). 

(iv) A discussion of how some assumptions ensure that the maximizing contin- 
uum does not touch the boundary of the underlying space except at a finite number 
of points. This ensures that variations of continua can be taken. 

(v) Proof of formula (22) involving the support of the equilibrium measure on 
the maximizing continuum and the external field (Theorem 5). 
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(vi) Proof that the support of the equilibrium measure on the maximizing con- 
tinuum consists of a union of finitely many analytic arcs. 

(v) Proof that the maximizing continuum is an S-curve (Theorems 7 and 8). 

The paper is organized as follows. In the rest of section 1, we introduce the 
appropriate topology for our set of continua that will provide the necessary com- 
pactness. In section 2, we prove a " Gauss- Frostman" type theorem which shows 
that the variational problem that we wish to solve is not vacuous. In section 3, we 
present the proof of upper semicontinuity of a particularly defined "energy func- 
tional" . In section 4, we present a proof of existence of a solution of the variational 
problem. Existence is thus derived from the semicontinuity and the compactness 
results acquired earlier. In section 5, we show that, at least under a simplifying 
assumption, the "max-min" solution of the variational problem does not touch the 
boundary of the underlying domain, except possibly at some special points. This 
enables us to eventually take variations and show that the max-min property im- 
plies regularity of the support of the solution and the S-property in sections 6 and 
7. By regularity, we mean that the support of the maximizing measure is a finite 
union of analytic arcs. In section 8, we conclude by stating the consequence of 
the above results in regard to the semiclassical limit of the nonlinear Schodingcr 
equation. 

We also include three appendices. The first one discusses in detail some topo- 
logical facts regarding the set of closed subsets of a compact space, equipped with 
the so-called Hausdorff distance. The fact that such a space is compact is vital 
for proving existence of a solution for the variational problem. The second appen- 
dix presents the semiclassical asymptotics for the initial value problem (9)-(10) in 
terms of theta functions, under the S-curve assumption (as in [3]). It is included 
so that the connection with the original motivating problem of semiclassical NLS 
is made more explicit. The third appendix shows how to get rid of the simplifying 
assumption introduced in section 5. 
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Following [6] (see Appendix A.l) we introduce an appropriate topology on F. 
We think of the closed upper half-plane H as a compact space in the Riemann 
sphere. We thus choose to equip H with the "chordal" distance, denoted by do, 
that is the distance between the images of z and Q under the stereographic pro- 
jection. This induces naturally a distance in K (so (io(0+,0_) ^ 0). We also de- 
note by do the induced distance between compact sets E,F in IK: do{E,F) = 
max z< zEminQ e Fdo{z, (). Then, we define the so-called Hausdorff metric on the set 
/(K) of closed non-empty subsets of IK as follows. 

(11) d K (A, B) = sup(do(A, B), d (B, A)). 

In appendix Al, we prove the following. 

LEMMA A.l. The Hausdorff metric defined by (11) is indeed a metric. The set 
/(K) is compact and complete. 

Now, it is easy to see that F is a closed subset of J(K). Hence F is also compact 
and complete. 

REMARKS. 

1. Because of the particular symmetry ip(x) = ip(—x) of the solution to the 
Cauchy problem (9)-(10) we will restrict ourselves to the case x > from now on. 
We then set J = 1 and the external field is 

(3a) 4>{z) = - J ' G(z;ri)dii°(ri) - Re[in J p°(v)dr] + 2i(zx + z 2 t)]. 

2. The function p° expresses the density of eigenvalues of the Lax operator 
associated to (9), in the limit as h — > 0. WKB theory can be used to derive an 
expression for p° in terms of the initial data tp°(x) via an Abel transform (see [3]), 
from which it follows that 

Re[p°(z)} = 0, for z e [0,iA], 
Im[p°(z)} > 0, for z G (0,iA]. 

The rest of the conditions (1) are not a necessary consequence of WKB theory. 
In particular, it is not a priori clear what the analyticity properties of p° are. In 
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this paper, we assume, for simplicity, that p° admits a continuous extension in 
the closed upper complex plane which is holomorphic in the open upper complex 
plane. We also assume that imp is positive in the real axis. This will be used 
later to show that the maximizing continuum docs not touch the real line, except at 
0+,0_,oo. It is a simplifying but not essential assumption. All conditions (1) are 
satisfied in the simple case where the initial data are given by ip(x,0) = Asechx, 
where A is a positive constant. 

3. It follows that is a subharmonic function in H which is actually harmonic 
in K ; it also follows that it is upper semicontinuous in H. It is then subharmonic 
and upper semicontinuous in H except at infinity. 

4. Even though in the end we wish that the maximum of E,p(X F ) over "continua" 
F is a regular curve, we will begin by studying the variational problem over the 
set of continua F and only later (in section 6) we will show that the maximizing 
continuum is in fact a nice curve. The reason is that the set F is compact, so once 
we prove in section 3 the upper semicontinuity of the energy functional, existence 
of a maximizing continuum will follow immediately 
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2. A GAUSS-FROSTMAN THEOREM 

We claim that for any continuum F e F, not containing the point oo and ap- 
proaching + , 0_ non-tangcntially to the real line, the weighted energy is bounded 
below and the equilibrium measure A F exists. This is not true for any external 
field, but it is true for the field given by (3a) because of the particular behavior of 
the function p° near zero. 

We begin by considering the equilibrium measure on the particular contour Fq 
that wraps itself around the straight line segment [0,iA], say X F . We have 

PROPOSITION 1. Consider the contour Fq e F consisting of the straight line 
segments joining 0+ to iA + = iA and iA = iA- to 0_. The equilibrium measure 
Xq exists. Its support is the imaginary segment [0,ibo(x)], for some < bo(x) < A, 
lying on the right of the slit [0,iA]. It can be written as p(z)dz where p(z) is a 
differentiable function in [0,ib(x)], such that p — p° belongs in the Holder class with 
exponent 1/2. 

PROOF: See section 6.2.1 of [3]; p(z) can be expressed explicitly. The Holder 
condition follows from (5.14) of [3]. Note that the field <f) is independent of time on 
F , so Xq is also independent of time. 

From Proposition 1, it follows that the maximum equilibrium energy over con- 
tinua is bounded below. 

(12) max Fe¥ E ljj {X F ) = max Fer min^ eM (F)E < j > (p) > -oo. 

The following formula is easy to verify. 



(13) 



(m) - E^X F ) = E(p — X F ) + 2 / (V xF + <j>)d{» - A F ), 



for any p which is a positive measure on the continuuum F. Here 




where again G(u, v) is the Green function for the upper half-plane. 
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To show that E$ (//) is bounded below, all we need to show is that the difference 
E^(fj,) — E t j > {X F ) is bounded below. 

Note that since V x " +0 = 0, on supp(X F ), the integral in (13) can be written 

as / {v xF + <f>)dfi. 

We have 

p rba(x) \ z 4-hi\ 

V x (z) + 4>(z)= / W-—±(-p) t=0 du + <t ) = 
Jo \z-iu\ 

= -Re[ [ b0(X) log^^lu^du] + 0(\z\) = 
Jo \z-iu\ 

= 0(\z\) near z = 0. 

So we can write V + <f> > c(A,x)\z\ in a neighborhood of z — 0, where c(A, x) 
will be some negative constant independent of z. Note that the dependence on t is 
not suppressed, but it is of order 0(|,z 2 |). 

It is now not hard to see that the 0(|^|) decay implies our result. 

Write fi — Ma, where M > is the total mass of n and a is a probability 
measure (on F). Choose e such that for \u\ < e we have V + 4> > c{A,x)\u\. 
Then 



E+in) - E^y ) > / G{u, v)d(u - A )(u)d(n - A )(«) 

( 14 ) r r 

+2/ (V A + (t>)(v)du(v) + 2c{A,x)\v\d{n-X F ){v). 

J\v\>t J\v\<e. 

The first integral of the right hand side (RHS) can be written as f\ u \> ei \ v \> e + 2 f\ u \ <e ,\ v \> e + I\u\,\v\<e ■ 
The sum of the first integral plus the second term of the RHS of (14) is bounded 
below, by the standard Gauss-Frostman theorem ([6], p. 135). It remains to consider 



(15) 

(/+/ )[/ G(u,v)d(n-X F )(u)]d(n-X F )(v)+ f 2c(A,x)\v\d(n-X F )(v) 

Jf J\v\>£ J\ll\<£ J\v\<e 



> 



([ + f )[[ G(u,v)d(Ma-X F )(u) + 2c(A,x)\v\]d(Ma-X F )(v). 

Jf J\v\>e J\u\<£ 



l\v\>e J\u\<t 

Now, it is easy to see that for M large, and since F is non-tangential to the real 
line, 



(16) / G{u 1 v)d{M(T- X t ){u) + 2c{A,x)\v\ > \vlogv\ + 2c(A,x)\v\. 

J\u\<e 



12 

By adjusting e if necessary, we have \vlogv\ > 2\c(A, x)\\v\, and hence the integral 
in (16) is positive. Integrating again with respect to Aider — d\ F , again for M large, 
we see that the integral of (15) is positive. 

Since when M is bounded above we have our estimates trivially, we get bound- 
cdness below over all positive M. 

We have thus proved one part of our (generalised) Gauss-Frostman Theorem. 

THEOREM 1. Let 4> be given by (3a). Let F be a continuum in K \ oo and 
suppose that F is contained in some sector it < a < arg(X) < f3 < as A — > 0. 
Let M(F) be the set of measures 6 M which are supported in F. (So, in particular 
their free energy is finite and <j> G Li(/z).) We have 

(17) inf l _ ie M(F)E <t> {n) > -oo. 

Furthermore the equilibrium measure on F exists, that is there is a measure X F e 
M(F) such that E4F] = £ (A F ) = inf^ M{F) E^{fi). 

PROOF: The proof that (17) implies the existence of an equilibrium measure is 
a well known theorem. For our particular field (f> given by (3) it is easy to prove. 
Indeed, the identity 

E{jjl -v) = 2E^(n) + 2E <p { V ) - 4£ (^) 

implies that any sequence fi n minimizing E^lp) is a Cauchy sequence in (un- 
weighted) energy. Since the space of positive measures is complete (see for example 
[7], Theorem 1.18, p. 90), there is a measure ^io such that E(fi n — fio) — > 0. We then 
have E(fi n ) — > E(fio) < +oo and hence fi n — > weakly (sec e.g. [7], p. 82-88; this 
is a standard result). 

The fact that <j> e £i(/io) is trivial for our particular field. 
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3. SEMICONTINUITY OF THE ENERGY FUNCTIONAL 



We consider the functional that takes a given continuum F to the equilibrium 
energy on this continuum: 



(18) E : F - E^F] = E^{\ F ) = inf (E{n) + 2 / 



and we want to show that it is continuous, if tp is continuous in ft. Note that this is 
not the case for the field </> given by (3a), since it has a singularity at oo; that field 
is only upper semicontinuous. We will see how to circumvent this difficulty later. 
For the moment, ip is simply assumed to be a continuous function in EL 

THEOREM 2. If ip is a continuous function in H\ oo then the energy functional 
defined by (18) is continuous at any given continuum F not containing the point 
oo. 

PROOF: Suppose G G F, also not containing oo, with d&(F, G) < d, a small 
positive constant. Let A = A^ be the equilibrium measure on F and = A^ be the 
equilibrium measure on G. 

We consider the Green's balayage of on F, say fi. Then supp fi G F and 

= V on F, 
J udfi = J udp 1 

for any function u that is harmonic in H \ F and continuous in H. 
Similarly consider A, the balayage of A to G. We trivially have 

E^[G] < Et(\), 
E^[F] <E^(fi). 

LEMMA 1. Suppose Q G F, p some positive measure supported in K and fi is 
the Green's balayage to Q. Then 

E(ji) = E(ji) - E Q c(n), 
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where Eqc(^) is the unweighted Green energy with respect to Q c = K \ Q. In 
particular, since unweighted energies are nonnegative, 



PROOF: The first identity follows from the fact that — vanishes on Q 
and the real line, and is harmonic in Q c and superharmonic in K. 

Integrating E(fi) = J V^djl = J Vdfi - J V^dfi = /(V" - V£ c )d/x = E(n) - 
Eqc(h). The proof of the Lemma follows. 

So, let u$ be a function harmonic in H \ F such that = ip on F and u$ = 
on <9H \ F. By the definition of balayage one has J ipdjl — J u^dfi. 

We have E,^ [F] < E^ (fi) = E(fi) + 2 J ^dfi < E([i) + 2 J i\)d\i + 2j{u^- ij))dn = 
E^n) + 2 f(u^ - ipW = E^[G] + 2 J (u^, - i>)dfx. 

In a small neighbourhood of F, Fd — {z : d(z, F) < d}, we have 



We assumed here that the equilibrium measures on continua near F are bounded 
above. This is easy to see. Suppose, first, that the point oo is not in F. Indeed, on 
the support of the equilibrium measure A, we have 



If the equilibrium measures on continua near F were unbounded, then so would be 
the potentials V x . But ip is definitely bounded near F. This contradicts the above 
equality. 

Now given y e Fd, choose z <G F such that \z — y\ = d. The above expression 
(19) is less or equal than 



CmaXyzpJu^iy) - ip(y) - u^(z) + ip{z)\ < o(l) + Cmax ye pju^(y) -u^(z)\. 
It remains to bound \u^(y) — u^{z)\ by an o(l) quantity. 



E(fi) < E(n). 



(19) 




+ V = o. 
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The next Lemma is due to Milloux and can be found in [8] . 

LEMMA 2. Suppose D is an open disc of radius R, with center zo; let y be a 
point in D, F a continuum in C, containing z , and be the connected component 
of D \ F containing y. Let w(z) be a function harmonic in fi such that 

w(z) = 0, z G Fndn, 

w(z) = l, z g 

Then w(y) < C{^^fl 2 . 
PROOF: See [8], p.347. 

Now, select a disc of radius d 1 / 2 , centered on z. We have \u^(y) — u${z)\ = o(l) 
on the part of F lying in the disc, while \u^(y)—u^(z)\ is bounded by some positive 
constant M on the disc boundary. 

LEMMA 3. Let fibea domain, dfl = F 1 U F 2 and 

wi = 0, z G Fi, 
= 1, F 2 ; 

w 2 = l, z e Fi, 
= 0, z e F 2 . 

Suppose m is harmonic in and 

u(z) < e, z G Fi, 
u(z) < M, z G F 2 . 

Then u(z) < ew 2 (z) + Mwi(z). 
PROOF: Maximum principle. 

Now, using Milloux's Lemma, we get \u^{y) — u^,(z)\ < o(i)w 2 (z) + Mw\(z) < 

1 /2 

o(l) + MC 1 -^ < o(l) + MCd 1 ' 2 . This concludes the proof of Theorem 2. 
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We now recall that the energy continuity proof was based on the continuity of 
ip. In our case, <j> is upper semicontinuous and discontinuous at oo. Still we can 
prove that the energy is upper semicontinuous and that will be enough. 

THEOREM 3. For the external field given by (3a), the energy functional defined 
in (18) is upper semicontinuous. 

PROOF: We first note that if the external field <fi' is upper semicontinuous away 
from infinity then so is the energy functional that takes a given continuum F to 
the equilibrium energy of F. Indeed, if <f>' is upper semicontinuous away from 
infinity, then there exists a sequence of continuous functions (away from infinity) 
such that <p n I <j)' . Each functional E^, n [F] is continuous, away from infinity, and 
E<p n [F] I Ep[F]. So, E^i[F\ is upper semicontinuous, away from infinity. 

Now consider the field <j) given by (3a). Let F be a continuum. If oo is not in F, 
then we're done. If oo £ F, let A = A F be the equilibrium measure. We can assume 
that on the equilibrium measure <j> is bounded by 0. Indeed, on the support of the 
equilibrium measure A, we have 

V x + </> = 0. 

But V x > 0, so 4> < 0. 

This means that we can change (f> to <f>' = min((f>, 0), which is an upper semicon- 
tinuous function. Theorem 3 is proved. 

REMARK. If we naively consider the functional taking a measure to its weighted 
energy we will see that it is not continuous even if the external field is continuous. 
It is essential that the energy functional is defined on equilibrium measures. 
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4. PROOF OF EXISTENCE OF A MAXIMIZING CONTINUUM 

THEOREM 4. For the external field given by (3a), there exists a continuum 
FeF such that the equilibrium measure X F exists and 

E <j> [F}(= E <j> (X F )) = max F€ wmin li€ M(F)E < i ) (fj,). 

PROOF: We know (see for example section 2) that there is at least one continuum 
F for which the equilibrium measure exists and E ( f,(X F ) > — oo, for all time. On 
the other hand, clearly E ( p(X F ) < for any F. Hence the supremum over continua 
in F is finite (and trivially nonpositive), since F is compact. Call it L. 

We can now take a sequence F n such that E<f,[F n ] — > L. Choose a convergent 
subsequence of continua F n — * F, say. By upper semicontinuity of the weighted 
energy functional, 

HmsupE4F n } <E 4> [F]<L = UmE^}. 
So L = E<f,[F\. The theorem is proved. 

5. ACCEPTABILITY OF THE CONTINUUM 

We have thus shown that a solution of the maximum-minimum problem exists. 
We do not know yet that the maximizing continuum is a contour. Clearly the 
pieces of the continuum lying in the region where the external field is positive do 
not support the equilibrium measure and by the continuity of the external field 
they can be perturbed to a finite union of analytic arcs. The real problem is to 
show that the support of the equilibrium measure is a finite union of analytic arcs. 
This will follow from the analyticity properties of the external field. 

Note that the maximizing continuum cannot be unique, since the subset where 
the equlibrium measure is zero can be perturbed without changing the energy. A 
more interesting question is whether the support of the equilibrium measure of the 
maximizing contour is unique. We do not know the answer to this question but it 
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is not important as far as the application to the semiclassical limit of the nonlinear 
Schrodinger equation is concerned. (See Appendix A2.) 

It is important however, that the maximizing continuum does not touch the 
boundary of the space IK except of course at the points 0~, + , and perhaps at oo. 
This is to guarantee that variations with respect to the maximizing contour can be 
properly taken. 

The proof of the acceptability of the continuum requires two things. 

(i) The continuum does not touch the real negative axis. 

(ii) The continuum does not touch the real positive axis. 
We will also make the folowing assumption. 

ASSUMPTION (A). The continuum maximizing the equilibrium energy does 
not touch the linear segment (0, 

REMARK. Assumption (A) is not satisfied at t — 0, where in fact the continuum 
is a contour Fq wrapping around the linear segment [0,iA]. However, the case 
t = is well understood. The equilibrium measure for Fq exists and its support 
is connected. On the other hand assumption (A) is satisfied for small t > 0. (See 
Chapter 6 of [3].) 

REMARK. It is conceivable that at some positive to there is an x for which 
assumption (A) is not satisfied. It can in fact be dropped but the analysis of the 
semiclassical limit of NLS will get more tedious; see Appendix A3. 

PROPOSITION 2. The continuum maximizing the equilibrium energy does not 
touch the real axis except at the points zero and possibly infinity. 
PROOF: (i) If z < 0, then 4>(z) = it /° \p°(v)\dr] > 0. 

This follows from an easy calculation, using the conditions defining p°. But 
we can always delete the strictly positive measure lying in a region where the 
field is positive and make the energy smaller. So even the solution of the "inner" 
minimizing problem must lie away from the real negative axis. 

(ii) If z > 0, then again a short calculation shows that d f^ z > 0, for t > 0. 

It is crucial here that if u € K then G(u, v) = 0, while if both u, v are off the real 
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line G(u, v) > 0. Hence, for any configuration that involves a continuum including 
points on the real line, we can find a configuration with no points on the real line, 
by pushing measures up away from the real axis, which has greater (unweighted 
and weighted) energy. So, suppose the maximizing continuum touches the axis. 
We can always push the measures up away from the real axis and end up with a 
continuum that has greater minimal energy, thus arriving at a contradiction. 
The proposition is now proved. 

REMARK. It is also important to consider the point at infinity. We cannot prove 
that the continuum does not hit this point. (In fact, our numerics ([3], Chapter 6) 
show that it may well do so.) In connection with the semiclassical problem (9)-(10) 
as analyzed in [3] , it might seem at first that the maximizing continuum should not 
pass through infinity. Indeed, the transformations (2.17) and (4.1) of [3] implicitly 
assume that the continuum C lies in C. Otherwise, one would lose the appropriate 
normalization for M at infinity. However, one must simply notice that infinity is 
just an arbitrary choice of normalizing point, once we view our Riemann-Hilbcrt 
problems in the compact Riemann Sphere. The important observation is that the 
composition of transformations (2.17) and (4.1) (which are purely formal, i.e. no 
estimates are required and no approximation is needed) docs not introduce any 
bad (essential) singularities. In the end, the asymptotic behavior of N" is still the 
identity as z — > oo in the lower half-plane and non-singular as z — > oo in the upper 
half-plane. So, in the end it is acceptable for a continuum to go through the point 
infinity. 
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6. TAKING SMALL VARIATIONS 

We now complexify the external field and extend it to a function in the whole 
complex plane, by turning a Green's potential to a logarithmic potential. We will 
thus be able to make direct use of the results of [5] . 

We let, for any complex z, 

/iA c-iA 
log{z - t])p°{r])dr] - (2ixz + 2itz 2 + in A 7 ?)*?) 
-iA Jz 

and Vr — ReV be the real part of V. In the lower half-plane the function p° is 
extended simply by 

(20a) pV) = (P° (I))'- 

Note right away that the field <f> defined in (3a) is the restriction of Vr — ReV to 
the closed upper half-plane. 

The actual contour of the logarithmic integral is chosen to be the linear segment 
X joining the points — iA, 0, iA. The branch of the logarithm function log(z — rj) is 
then defined to agree with the principal branch as z — > oo , and with jump across 
the very contour \. 

The unweighted Green's energy (4) can be written as 

(21) E Vr {h)= f [ l og —}—d i i{u)dp{v) + 2[ V R (u)dp(u), 
where the measures \x are extended to the lower half complex plane by 

(21a) mCO = -/*(*)■ 

(So they are "signed" measures.) 

Having established in section 5 that the contour solving the variational problem 
does not touch the real line we can take small variations of measures and contours, 
never intersecting the real line, and keeping the points 0_,0 + fixed. In view of 
(21a) we can think of them as variations of measures symmetric under (21a) in the 
full complex plane, never intersecting the real line, and keeping the points 0_,0 + 
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fixed. The perturbed measures do not change sign. The fact that oo can belong to 
the contour is not a problem. Our variations will keep it automatically fixed. 

The first step is to show that the solution of the variational problem satisfies a 
crucial relation. 

REMARK. It is not hard to see that the variational problem of Theorem 4 is 
actually equivalent to the variational problem of maximizing equilibrium measures 
on continua in the whole complex plane, under the symmetry (21a) and the condi- 
tion that measures are to positive in the upper half-plane and negative in the lower 
half-plane. 

THEOREM 5. Let F be the maximizing continuum of Theorem 4 and X F be 
the equilibrium measure minimizing the weighted logarithmic energy (6) under the 
external field Vr — ReV where V is given by (20). Let fi be the extension of \ F to 
the lower complex plane via (i(z*) = — u(z). Then 



PROOF: We first need to prove the following. 

THEOREM 6. Let L be a critical point of the functional taking a continuum 
L e F to Ev n (^ r ), and assume that F is not tangent to R. Also assume that T 
does not touch the segment [0, iA] except at zero. Let \x be the extension of A r via 
n(z*) — —u(z), Or be an open set containing the interior of Tur* and h G C 1 (Or) 
such that h(0) = 0. We have 



PROOF: Consider the family of (signed) measures {/i T ,r G C, \t\ < tq} defined 
by d/j, T (z T ) — du(z) where z T = z + rh(z)), or equivalently, f f(z)du T (z) = 
J f(z T )dfi(z), f G Li(Or)- Assume that r is small enough so that the support 
of the deformed continuum does not hit the linear segment (0, iA] or a non-zero 
point in the real line. 



(22) 




(23) 
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With h = h(u,v) = h{u) u Z h v (v \ we have = 1 + rh, so that log^^ - 

l °9^h\ = + rh\ = -Re(rh) + 0(r 2 ). 

Integrating with respect to dp,(u)dfi(v) we arrive at 

(24) EipJ) - E(fi) = —Re [r J hdp,(u)dfj,{v)} + 0(t 2 ), 

where E{p) denotes the free logarithmic energy of the measure p,. Also, 
J V R dfi T - J V R dii = 2 J (V R (u T ) - V R {u))d^{u) = 
2Re[r J V ' (u)h(u)d^u)} + 0(t 2 ). 

Combining with the above, 

(25) E VR (fi T ) - E Vr (/j,) = Re(—T J hd^{u)d^{v) + 2t J V'hd^) + 0(t 2 ). 

So, if fi is (the symmetric extension of) a critical point of the map fi — > Ey R {p) 
the linear part of the increment is zero. In other words given a C 1 function h and 
a measure p, the function Ey R {p T ) of r is differentiable at t = and the derivative 
is 

(26) Re(-H(p)), where H(p) = j j hdp 2 -2 J V'hd^. 

But what we really want is the derivative of the energy as a function of the equi- 
librium measure. This function can be shown to be differentiable and its derivative 
can be set to zero at a critical continuum. 

Indeed, we need to show the following. 

LEMMA 4. 

±E VR ((X r r)\r=o = ^E Vr (\ t t)\t=o = 0. 

In the relation above T T = supp(X r ) T . The first derivative is of a function of general 
measures. The second derivative is of a function of equilibrium measures. 

PROOF: Define the measure er T with support T and such that (cr T ) T = X Tt . 

LEMMA 5. With H defined by (26), we have 

ReH(a T ) -► ReH(X r ), 
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as T — > 0. 

PROOF. By (25)-(26), we have 

E VR ((\ r r) E Vr (X t ) = -Re(rH(X T ) + 0(r 2 )), 

E Vn {\ T -) - E v Ja T ) = -Re(rH(a T ) + 0(r 2 )). 
On the other hand, E Vr {<j t ) > E Va (X r ), and E Vn {{X r ) T ) > E Vn (X r -). It follows 

that 

E VR (a T ) - Re(rH(a T )) + 0(r 2 ) = E Vr (X Tt ) < E Vr ((X t Y) = E VR (X r ) - i?e(rff(A r )) + 0(r 2 ). 

Hence E Vr (ct t ) -> E Va (X r ). 

As in the proof of Theorem 1, it follows that <r T — > A r weakly; see [7], pp. 82-88. 
It then follows immediately that H(a T ) — > H(X r ). This proves Lemma 5. 

To complete the proof of Lemma 4, we note that > Ev R ((X r ) T ) — Ey R {X T ) > 
e v r (>F t )-E Vr (<t t ) = -Re(TH(<r T )) + 0(T 2 ). Hence the derivative of E Va ((X r ) T ) 
at r = is equal to the derivative of Ev r {X Tt ) at r = which is equal to ReH(X r ). 
This proves Lemma 4 and Theorem 6. 

PROOF OF THEOREM 5. Consider the Schiffer variation, i.e. take h(u) = 

' V / u—z 

where z is some fixed point not in T. Note that h(0) = so that the deformation 
z T = z + rh(z) keeps the points 0+, 0- fixed. Also assume that r is small enough so 
that the support of the deformed continuum does not hit the linear segment (0, iA] 
or a non-zero point in the real line. We have 

h = h(u,v)= h ^- h ^=l- 



u — v (u — z)(v — z)' 

and therefore 

dfi(u) 2 



Next 



Re[ / / h(u, v)dfi(u)dfi(v)] = Re( / dfi(u)d/j,(v) — z 2 [ 

J J J J J su 

Re[ j 2V'(u)h(u)dn(u)} 
= 2Re[ J(u + z)V\uW{u) + z 2 J Y'to*^) , 



u — z 



= Rz[j 2V'(u)(u + z)dfj,(u)+ 

2z2 [v»-v'(*) Mu) + 2z , v , {z) /W], 

J u— z J u— z 



Theorem 5 now follows from Theorem 6. 
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REMARK. If our continuum is allowed to touch the point iA (so we slightly 
weaken assumption (A)) then we may need to keep points ±iA fixed under a small 

2 / 2. a 2 \ 

variation. We can then choose the Schiffer variation h(u) — - u _ z — -. We will 
arrive at a similar and equally useful formula. 

In general if one wants to keep points a±, ...,a s fixed, the appropriate Schiffer 
variation is h(u) = IIi =^"~ a '' > . 

PROPOSITION 3. The support of the equilibrium measure consists of a finite 
number of analytic arcs. 

PROOF: Theorem 5 above implies that the support of [i is the level set of the 
real part of a function that is analytic except at countably many branch points. In 
fact, supp^i is characterized by / log \ u \ z \ d/i(u) + Vr(z) = . From Theorem 5 we 
get 

(27) Re[J^ + V(z)]=Be[(R^(z)m 



where 



(28) 



suppfi z u 



+\([ 2(u + z)V'(u) dfi(u)). 

Z J suppp. 



This is a function analytic in K, with possibly a pole at zero. By integrating, we 
have that supp(p) is characterized by 

(29) Re J (Rf,) 1/2 dz = 0. 

The locus defined by (29) is a union of arcs with endpoints at zeros of To 
see this, consider (29) as an equation of two real variables f(u,v) = and try to 
solve for u as an analytic expression of v. One can only do this if the derivative 
of / is nonzero, which means (via the Cauchy-Riemann equations) that R^ ^ 0. 
The points where i? M = are exactly the points where the analytic arcs cannot be 
continued. 
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By conditions (1) for p°, R^ is blowing up at the point oo (at least for t > 0; 
but the case t = is well understood: the equilibrium measure consists of a single 
analytic arc; see section 5). Hence it can only have finitely many zeros near infinity, 
otherwise they would have to accumulate near oo and then would be there. On 
the other hand, R^ cannot have an accumulation point of zeros at z = 0, because 
even if the pole at were removed (the coefficients of \, - being zero), i? M would 
be holomorphically extended across z = 0. So, R^ can only have a finite number 
of zeros in K. It follows that the support of the maximizing equilibrium measure 
consists of only finitely many arcs. 

REMARK. Of course, conditions (1) can be weakened. We could allow p° to 
have a pole at infinity of order other than two. But our aim here is not to prove the 
most general theorem possible, but instead illustrate a method that can be applied 
in the most general settings under appropriate amendments. 

REMARK. The assumption that p° is continuous and hence bounded at infinity 
is only needed to prove the finiteness of the components of the support of the 
equilibrium measure of the maximizing continuum. If it is dropped then we may 
have an infinite number of components for isolated values of x, t. This will result 
in infinite genus representations of the semiclassical asymptotics. Of course infinite 
genus solutions of the focusing NLS equation are known and well understood. So 
the analysis of [3] is expected to also apply in that case, although it will be more 
tedious. 

For a justification of the "finite gap ansatz", concerning the semiclassical limit 
of focusing NLS, it only remains to verify the " S-property" . 



Note that 



(30) 
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7. THE S-PROPERTY 
THEOREM 7. (The S-property) 

Let C be the contour maximizing the equilibrium energy, for the field given by 
(3a) with conditions (1). Let fj, be the extension of its equilibrium measure to the 
full complex plane via (21a). Assume for simplicity that assumption (A) holds. Let 
= f mpPI >l°9(th)M«), X R {z) = ReX(z) = f suppfl log(^W(u),W„ = 
X'. Then, at any interior point of supp^i other than zero, 

(8a) ^ + *«) = ^ + ^)< 

where the two derivatives above denote the normal derivatives, on the + and — 
sides respectively. 

PROOF: From Theorem 5, we have 

\Re(W,(z) + V'(z))\ = \Re(R,)\ 1 / 2 . 

Using the definition for X, the above relation becomes 

\±Re(X + V)\ = \Re(R„)\ 1 / 2 . 

Now, Re(X+V) = on the support of the equilibrium measure. So, in particular 
Re(X + V) is constant along the equilibrium measure. Hence \£Re{X + V)\ must 
be equal to the modulus of each normal derivative across the equilibrium measure. 
So, 

\-^(V R + ReX)\ = \-^Re(X + V)\ = \Re{R,)\ 1/2 - 

Hence, 

\jL(V R + ReX)\ = \J-(V R + ReX)\. 

an + an- 

But it is easy to see that both LHS and RHS quantities inside the modulus sign 
are negative. This is because Vr + ReX = on suppfj, and negative on each side 
of suppfi. Hence result. 
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REMARK. Once Theorem 7 is proved it follows by the Cauchy-Riemann equa- 
tions that (Vf + ImX) + + (Vj + ImX)_ is constant on each connected component 
of suppn, which means that Imcp is constant on connected components of the con- 
tour, where <fi is as defined in formula (4.13) of [3]. This proves the existence of the 
appropriate " g- functions" in [3]. 

We recapitulate our results in the folllowing theorem, set in the upper complex 
half- plane. Note that (8a) is the "doubled up" version of (8). 

THEOREM 8. Let <f> be given by (3a), where p° satisfies conditions (1). Under 
assumption (A), there is a piecewise smooth contour C e F, containing points 
0+, 0- and otherwise lying in the cut upper half-plane K, with equilibrium measure 
X c , such that supp(X ) consists of a union of finitely many analytic arcs and 

E<p(X c ) = maxcewE^X ) = maxc'er[inf li eM(F)E^ ((/,)]. 
On each interior point of supp(\ c ) we have 

where is the Green's potential of the equilibrium measure X c (see (5)) and the 
two derivatives above are the normal derivatives. 

A curve satisfying (8) such that the support of its equilibrium measure consists 
of a union of finitely many analytic arcs is called an S-curve. 

PROOF: The fact that the maximizing continuum C is actually a contour is 
proved as follows. If this were not the case, then we could choose a subset of C, say 
F, which is a contour, starting at 0+ and ending at 0_, and going around the point 
iA. Clearly, by definition, the equilibrium energy of C is less than the equilibrium 
energy of F, i.e. E c f ) (X c ) < £ (A F ). On the other hand, since C maximizes the 
equilibrium energy, we have E^{\ F ) < E^,(X C ). So E^(X F ) — E^{\ c ). 
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8. CONCLUSION. 

In view of the interpretation of the variational problem in terms of the semiclas- 
sical NLS problem, we have the following result. 

Consider the semiclassical limit (h — ► 0) of the solution of (9)-(10) with bell- 
shaped initial data. Replace the initial data by the so-called soliton ensembles data 
(as introduced in [3]) defined by replacing the scattering data for ip(x,0) = ipn(x) 
by their WKB-approximation, so that the spectral density of eigenvalues is 

dpf KB (r ] ) := p°(v)X[0M](v)dv + pV)*X[-iA,o] 

with p\n) := ^ ; dX = -4- VM*) 2 + V 2 dx, 

I" Jx-( V ) V A(x) 2 +i] 2 ndr] J x _ (v) 

for rye (0, iA), where x_ (77) < x + (rf) are the two real turning points, i.e. (A(x±)) 2 + 
i] 2 = 0, the square root is positive and the imaginary segments (—iA, 0) and (0, iA) 
are both considered to be oriented from bottom to top to define the differential dr}. 

Assume that p° satisfies conditions (1). Then, under assumption (A), asymp- 
totically as h — > 0, the solution tp(x,t) admits a "finite genus description". (For a 
more precise explanation, see Appendix A2.) 

The proof of this is the main result of [3], assuming that the variational problem 
of section 1 has an S-curve as a solution. But this is now guaranteed by Theorem 
8. 

REMARK. For conditions weaker than the above, the particular spectral density 
p° arising in the semiclassical NLS problem can conceivably admit branch singular- 
ities in the upper complex plane and condition (1) will not be satisfied. We claim 
that even in such a case the finite gap genus can be justified, at least generically. 
The proof of this fact will require setting the variational problem on a Riemann 
surface with moduli at the branch singularities of p°. 



EXISTENCE AND REGULARITY FOR AN ENERGY MAXIMIZATION PROBLEM IN TWO DIMENSION89 

REMARK. Consider the semiclassical problem (9)- (10) in the case of initial data 
tpo(x) = Asechx, where A > 0. Then the WKB density is given by p° — i (sec (3.1) 
and (3.2) of [3]; note that condition (1) is satined). So the finite genus ansatz holds 
for any x, t, as long as the assumption (A) of section 5 holds. But then assumption 
(A) can be eventually dropped; see Appendix A3. 

REMARK. The behavior of a solution of (9) in general depends not only on 
the eigenvalues of the Lax operator, but also on the associated norming constants 
and the reflection coefficient. In the special case of the soliton ensembles data the 
norming constants alternate between —1 and 1 while the reflection coefficient is 
by definition zero. More generally, for real analytic data decaying at infinity the 
reflection coefficient is exponentially small everywhere except at zero and can be 
neglected (although the rigorous proof of this is not trivial). 
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APPENDIX Al. COMPACTNESS OF THE SET OF CONTINUA 

In this section we prove that the sets J(K) and hence F defined in section 1 are 
compact and complete. 

As stated in section 1, the space we must work with is the upper half-plane: 
H = {z : Imz > 0}. The closure of this space is H = {z : Imz > 0} U {oo}. Also 
K = {z : Imz > 0} \ {z : Rez — 0, < Imz < A}. In the closure of this space, K, 
we consider the points ix + and ix-, where < x < A as distinct. 

Even though we eventually wish to consider only smooth contours, we are forced 
to a priori work with general closed sets. The reason is that the set of contours is 
not compact in any reasonable way, so it seems impossible to prove any existence 
theorem for a variational problem defined only on contours. Instead, we define F 
to be the set of all "continua" F in K (i.e. connected compact sets, containing the 
points 0+, 0_). 

Furthermore, we need to introduce an appropriate topology on F, that will make 
it a compact set. In this we follow the discussion of Dicudonne ([6], chapter III. 16). 

We think of the closed upper half-plane H as a compact space in the Riemann 
sphere. We thus choose to equip H with the "chorda!" distance, denoted by do(z, (), 
that is the distance between the images of z and £ under the stereographic projec- 
tion. This induces naturally a distance in K (so, for example, do(0+,0_) =/= 0). 
We also denote by do the induced distance between compact sets E, F in K: 
do(E, F) — max ze Emin^ e FdQ(z, £). Then, we define the so-called Hausdorff metric 
on the set /(K) of closed nonempty subsets of K as follows. 

(A.l) d K (A, B) = sup(d Q (A, B), d (B, A)). 

LEMMA A.l. The Hausdorff metric defined by (A.l) is indeed a metric. The 
set /(K) is compact and complete. 

PROOF: It is clear that d^(A,B) is non-negative and symmetric by definition. 
Also if d^{A,B) = 0, then do(A, B) = 0, hence max ze Amin^ & Bdo(zX) — and 
thus for all z e A, we have min^ e Bdo{z, () = 0. In other words, z G B. By 
symmetry, A = B. 
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The triangle inequality follows from the triangle inequality for do- Indeed, sup- 
pose A,B,C G I(K). Then d K (A,B) = sup(d (A,B),d (B,A)) = d (A,B), with- 
out loss of generality. Now, 

d (A, B) = max zeA minQ eB d a (z, () 

< max zeA min c <= B min Coe c(do(z, Co) + d (Co, C)). 
by the triangle inequality for do. Let z = zq G A be the value of z that maximizes 
min ieB min (oe c(do(z, Co) +d (Co,C))- Tms is then 

min CeB min CoeC (do(z , Co) + d (Co, 0) 

< min CoeC do{z , Co) + min CeB min CoeC d {(oX) 
< max ze Amin (oeC do{z, Co) + max ieB min ioeC do(Co, C) 

< do(A, C) + d (B, C) < d K (A, C) + d K (B, C). 
The result follows from symmetry. 

We will next show that /(K) is complete and precompact. Since a prccompact, 
complete metric space is compact ([6], proposition (3.16.1)) the proof of Lemma 
A.l follows. 

LEMMA A. 2 If the metric space E equipped with a distance do is complete, then 
so is /(E), the set of closed nonempty subsets of E , equipped with the Hausdorff 
distance 

d m (A, B) = sup(d (A, B), d Q (B, A)), 

for any closed nonempty subsets A, B, where do(A, B) = max a& AminbeBdo(a, b). 
Furthermore, if E is precompact, then so is /(E). 

PROOF: Suppose E is complete. Let X n be a Cauchy sequence in /(E). We will 
show that X n converges to X = n rl >oU p >oA n+p . (Overbar denotes closure.) 
Indeed, given any e > 0, 

d (X n ,X) = max x( zx n rnin ve xdn(x,y) 

< max xe x n max ye u p > x n+p do(x,y) < e, 
for large n, by the completeness of E. Similarly, 

d (X,X n ) = max x&x min ye x n do(x,y) 

< max xe u p > x n+p rnin ye x n do( x ,y) < e - 
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Next, suppose E is precompact. Then, by definition, given any e > 0, there is a 

finite set, say S — {si, S2, , s n }, where n is a finite integer, such that any point 

x of E is at a distance d less than e to the set S. Now, consider the set of subsets 
of S, which is of course finite. Clearly every closed set is at a distance less than e 
to a member of that set: 

do(A, S) = max aeA min seS do(a, s) < e, 

d (S, A) = max se smin aeA d (a, s) < e, 

for any closed nonempty set A. Hence dE{A, S) < e. 

So, any closed nonempty set A is at a distance less than e to the finite power set 
of 5\ So /(E) is precompact. 
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APPENDIX A2. THE DESCRIPTION OF THE SEMICLASSICAL LIMIT OF 
THE FOCUSING NLS EQUATION UNDER THE FINITE GENUS ANSATZ 

We present one of the main results of [3] on the semiclassical asymptotics for 
problem (9)-(10), in view of the fact that the finite genus ansatz holds. In partic- 
ular, we fix x, t and use the result that the support of the maximizing measure of 
Theorems 4 and 8 consists of a finite union of analytic arcs. 

First, we define the so-called g- function. Let C be the maximizing contour of 
Theorem 4. A priori we seek a function satisfying 



g{X) is independent of h. 



g{\) is analytic forX G C \ (C U C*). 



g(X) — > as X — > oo. 



g (A) assumes continuous boundary values from both sides of C U C* , 



denoted by g + {gJ) on the left (right) of C U C* . 



g(X*) + g(X)* = for all X e C \ (C U C*). 



The assumptions above are satisfied if we write g in terms of the maximizing 



equilibrium measure of Theorem 8, dp — dX c = p(n)dn, doubled up according to 



(21a). Indeed, 




for an appropriate definition of the logarithm branch (see [3]). 



For A G C, define the functions 



9(X) :=i(g + (X)-g_(X)), 




iA 








■iA 



where p (rf) is the holomorphic function (WKB density of eigenvalues) introduced 
in section 1 (see conditions (1)). 



The finite genus ansatz implies that for each x, t there is a finite positive integer 



G such that the contour C can be divided into "bands" [the support of p(n)dn] and 
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"gaps" (where p = 0). We denote these bands by Ij. More precisely, we define 
the analytic arcs Ij,I*,j = l,...,G/2 as follows (they come in conjugate pairs). 
Let the points Xj, j = 0, G, in the open upper half-plane be the branch points 
of the function g. All such points lie on the contour C and we order them as 
Ao, Ai, Ag, according to the direction given to C. The points A2,AJ,...,Aq are 
their complex conjugates. Then let Jo = [0, Ao] be the subarc of C joining points 
and A . Similarly, Ij — [\ 2 j-i, A 2 j], j = l,...,G/2. The connected components of 
the set C \ Lij(Ij U IJ) are the so-called "gaps", for example the gap Fi joins A to 
Ai, etc. 

It actually follows from the properties of g,p that the function 0(A) defined on 
C is constant on each of the gaps Tj, taking a value which we will denote by 6j, 
while the function $ is constant on each of the bands, taking the value denoted by 
aj on the band Ij. 

The finite genus ansatz for the given fixed x, t implies that the asymptotics of 
the solution of (9)- (10) as h — * can be given by the next theorem. 

THEOREM A.l. Let x 0l t be given. The solution ip{x,t) of (9)-(10) is asymp- 
totically described (locally) as a slowly modulated G + 1 phase wavetrain. Setting 
x = x + hi and t = t + hi, so that xo, to are "slow" variables while x, t are "fast" 
variables, there exist parameters 

a,U = (U ,Ui,....,U G ) T , k = (fco.fci, ,k G ) T , w = (w ,wi, ,wg) T , Y = 

(10,^1, , Yg) t , Z = (Z ,Zi, ,Zg) t depending on the slow variables xq 

and to (but not x, t) such that 

il>(x,t) = ip(x + tub, to + hi) - a{xo,to)e iUo(xoM/h e iiko(x ^ to)£ - Wo{x ^ to)i) 
( A - 2 ) Q(Y(x ,t ) + iU(x ,to)/h + i(k(x ,t )x - w(x ,t )t)) 

Q(Z(xq, h) + iU(x ,to)/h + i(k(xo,to)x - w(x ,to)i)) 
All parameters can be defined in terms of an underlying Riemann surface X. 

The moduli of X are given by Aj, j = 0,....,G and their complex conjugates 

A*, j = 0, G. The genus of X is G. The moduli of X vary slowly with x, t, i.e. 

they depend on x ,io but not x,t. For the exact formulae for the parameters as 

well as the definition of the theta functions we present the following construction. 
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The Riemann surface X is constructed by cutting two copies of the complex 
sphere along the slits I (J Iq , Ij , I* , j = 1...., G, and pasting the "top" copy to the 
" bottom" copy along these very slits. 

We define the homology cycles aj,bj, j = 1,...,G as follows. Cycle a\ goes 
around the slit Iq U 7q joining Ao to Aq, remaining on the top sheet, oriented coun- 
terclockwise, 02 goes through the slits 7_i and I\ starting from the top sheet, also 
oriented counterclockwise, az goes around the slits 7_i, Iq U7q j h remaining on the 
top sheet, oriented counterclockwise, etc. Cycle 61 goes through I and 7i oriented 
counterclockwise, cycle hi goes through 7_i and 7i, also oriented counterclockwise, 
cycle &3 goes through 7_i and 7 2 , and around the slits 7_i,7 U Iq,Ii, oriented 
counterclockwise, etc. 

On X there is a complex G-dimensional linear space of holomorphic differentials, 
with basis elements Vk{P) for k = 1, . . . , G that can be written in the form 

G-i 



"* {P) = S '°R X{ P) " A(P) - 
where Rx{P) is a "lifting" of the function i?(A) from the cut plane to X: if P is 
on the first sheet of X then Rx{P) = R(X(P)) and if P is on the second sheet of 
X then Rx(P) = —R(X(P)). The coefficients Cfcj are uniquely determined by the 
constraint that the differentials satisfy the normalization conditions: 

Vk{P) = 2iridjk- 

From the normalized differentials, one defines a G x G matrix H (the period matrix) 
by the formula: 

Hjk = f v k {P). 



It is a consequence of the standard theory of Riemann surfaces that H is a sym- 
metric matrix whose real part is negative definite. 
In particular, we can define the theta function 



Q(w) := ^ cxp(in T i7n + n T w), 
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where H is the period matrix associated to X. Since the real part of H is negative 
definite, the series converges. 

We arbitrarily fix a base point Pq on X. The Abel map A : X — > Jac(X) is then 
defined componentwise as follows: 



A k (P;Po) ■= f MP'), k = l,...,G, 

JPn 



where P' is an integration variable. 

A particularly important clement of the Jacobian is the Riemann constant vector 
K which is defined, modulo the lattice A, componentwise by 



P^kk 

K k := m + 



2 2m 

where the index k varies between 1 and G. 

Next, we will need to define a certain meromorphic differential on X . Let Vl(P) 
be holomorphic away from the points ooi and oo 2 , where it has the behavior 

n(p) = d P (\(p)) + ^±y p^oox, 

O(P) = -dp(X(P)) + O (j^) , P - oo 2 , 
and made unique by the normalization conditions 

n(p) = o,j = i,...,G. 



Here p is a polynomial, defined as follows. 

First, let us introduce the function R(X) defined by 

k=a 

choosing the particular branch that is cut along the bands l£ and 1^ and satisfies 

hm - 7 t—- = -1. 
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This defines a real function, i.e. one that satisfies R(X*) = R(X)* . At the bands, we 
have R+(X) = — i?_(A), while R(X) is analytic in the gaps. Next, let us introduce 
the function k(X) defined by 

Next let 

H(X) = k(X)R(X). 
The function k satisfies the jump relations 

fc + (A)-fc_(A) = -_*2_ Aer+ur- 

fc + (A) - fc_(A) = -^y, A e J+ U J", 

and is otherwise analytic. It blows up like (A — A„) -1 / 2 near each endpoint, has con- 
tinuous boundary values in between the endpoints, and vanishes like 1/A for large 
A. It is the only such solution of the jump relations. The factor of R(X) renor- 
malizes the singularities at the endpoints, so that, as desired, the boundary values 
of H(X) are bounded continuous functions. Near infinity, there is the asymptotic 
expansion: 

H(X) = H G X G + Hg^X - 1 + --- + H 1 X + H q + OiX- 1 ) 

(A.3) 

= p(X) + 0(X- 1 ), 

where all coefficients Hj of the polynomial p(X) can be found explicitly by expanding 
i?(A) and the Cauchy integral k(X) for large A. It is easy to sec from the reality of 
9j and aj that p(X) is a polynomial with real coefficients. 

Thus the polynomial p(X) is defined and hence the meromorphic differential fi(-P) 
is defined. 

Let the vector U G C G be defined componentwise by 

Uj := / O(P). 
J ^ 

Note that O(P) has no residues. 
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Let the vectors V\ , V 2 be defined componentwise by 

V hk = (A k (Xl + ) + A k (X 2+ ) + A k {X* 3+ ) + ■■■ + A k {X G+ )) + A k {oo) + m + 
V 2 . k = (A k (X* 1+ ) + A k (X 2+ ) + A k (X* + ) + ■■■ + A k (X G+ )) - A k (oo) + ™ + 

where k = 1, G, and the + index means that the integral for A is to be taken on 
the first sheet of X, with base point A+. 
Finally, let 

« = §if>~ 

k n = d x U n , w n = -d t U n , n = 0,...,G, 

where 

Y = -A(oo) - Vi, Z = A(oo) - V u 

and Uo = — {9\ + «o) where 0\ is the (constant in A) value of the function 9 in the 
gap Ti and a is the (constant) value of the function <f) in the band I . 

Now, the parameters appearing in formula (A. 2) are completely described. 

We simply note here that the Ui and hence the fcj and Wi are real. We also note 
that the denominator in (A. 2) never vanishes (for any xq, t , x, t). 

REMARK. Theorem A.l presents pointwise asymptotics in x, t. In [3], these are 
extended to uniform asymptotics in certain compact sets covering the x,t-plane. 
Error estimates are also given in [3]. 

REMARK. As mentioned above, we do not know if the support of the equilibrium 
measure of the maximizing continuum is unique. But the asymptotic formula (A. 2) 
depends only on the endpoints of the analytic subarcs of the support. Since the 
asymptotic expression (A. 2) must be unique, it is easy to see that the endpoints 
also must be unique. Different Riemann surfaces give different formulae (except of 
course in degenerate cases: a degenerate genus 2 surface can be a pinched genus 
surface and so on). 
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APPENDIX A3. DROPPING ASSUMPTION (A) OF SECTION 5. 

In section 5, we have assumed that the solution of the problem of the maximiza- 
tion of the equilibrium energy is a continuum, say F, which does not intersect the 
linear segment [0,iA] except of course at 0+,0_. We also prove that F does not 
touch the real line, except of course at and possibly oo. This enables us to take 
variations in section 6, keeping fixed a finite number of points, and thus arrive at 
the identity of Theorem 5, from which we derive the regularity of F and the fact 
that F is, after all, an S-curve. 

In general, it is conceivable that F intersects the linear segment [0, iA] at points 
other than + , 0_. If the set of such points is finite, there is no problem, since we 
can always consider variations keeping fixed a finite number of points, and arrive 
at the same result (see the remark after the proof of Theorem 5). 

If, on the other hand, this is not the case, we have a different kind of problem, 
because the function V introduced in section 6 (the complexification of the field) 
is not analytic across the segment [— iA, iA]. 

What is true, however, is that V is analytic in a Riemann surface consisting of 
infinitely many sheets, cut along the line segment [— iA, iA]. So, the appropriate, 
underlying space for the (doubled up) variational problem should now be a non- 
compact Riemann surface, say L. 

Compactness is crucial in the proof of a maximizing continuum. But we can 
compactify the Riemann surface L by compactifying the complex plane. Let the 
map C — > L be defined by 

y = log(z — iA) — log(z + iA). 

The point z = iA corresponds to infinitely many y-points, i.e. y = —oo + i0, 9 £ K, 
which will be identified. Similarly, the point z = —iA corresponds to infinitely 
many points y = +oo + id, e K, which will also be identified. The point G C 
corresponds to the points km, k odd. 

By compactifying the plane we then compactify the Riemann surface L. The dis- 
tance between two points in the Riemann surface L is defined to be the correspond- 
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ing stereographic distance between the images of these points in the compactificd 
C. 

With these changes, the proof of the existence of the maximizing continuum in 
sections 1, 3, 4 goes through virtually unaltered. In section 6, we would have to 
consider the complex field V as a function defined in the Riemann surface L and 
all proofs go through. The corresponding result of section 7 will give us an S-curve 
C in the Riemann surface L. We then have the following facts. 

Consider the image D of the closed upper half-plane under 

y = log(z — iA) — log{z + iA). 

Consider continua in D containing the points y — ni and y = —ni. Define the 
Green's potential and Green's energy of a Borel measure by (4), (5), (6) and the 
equilibrium measure by (7). Then there exists a continuum C maximizing the 
equilibrium energy, for the field given by (3) with conditions (1). C does not touch 
<9B except at a finite number of points. By taking variations as in section 6, one 
sees that C is an S-curve. In particular, the support of the equilibrium measure on 
C is a union of analytic arcs and at any interior point of supp^i 

where the two derivatives above denote the normal derivatives. 

As far as the consequences of the above facts regarding the semiclassical limit 
of NLS, some more work on the details is necessary. Indeed, the analysis of the 
dynamics of the semiclassical problem (9)-(10) presented in [3] assumes that the 
S-curve C lies entirely in K U {0 + , 0_}. However, the explicit computation of the 
equilibrium measure on the S-curve and the derivation of the equations implicitly 
defining the endpoints of the components of the support of the equilibrium measure 
also make sense when extended to the Riemann surface L. As a result the statement 
of Theorem A.l is correct, when interpreted correctly, i.e. allowing for all contours 
involved in the definition of the line integrals appearing in formula (A. 2) as well as 
the bands Jo U Iq, Ij, Ij ,j = 1--,G and the gaps Tj to lie in L. 
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We plan to describe the details of the evolution of the S-curve in the Riemann 
surface L in a later publication. In particular, one must check that all the defor- 
mations described in [3] are valid; this is indeed true, as one can see by the analysis 
of [3] transferred to the Riemann surface L. 

REMARK. The "discrete-to-continuous" passage. 

The following caveat has to be addressed here. One of the problems dealt with 
in [3] was to transform the so-called " discrete" Riemann-Hilbert problem (problem 
4.1.1 of [3]) to a "continuous" Riemann-Hilbert problem (problem 4.2.1 of [3]). The 
initial (discrete) problem posed is to construct a (matrix) meromorphic function 
from the information on its poles. These poles indeed accumulate (as h goes to 
0) along the segment (-iA, iA). Then one transforms the problem to a properly 
speaking (continuous) holomorphic Riemann-Hilbert problem with jumps, by con- 
structing two contours, one encircling the poles with positive imaginary part and 
the other encircling the poles with negative imaginary part, and redefining the ma- 
trix functions inside the contours so that the singularities are removed. Naturally 
the information concerning the poles appears in the arising jump along the two 
loops. 

One then (formally at first) approximates the jump on the loops in an obvious 
way (it is a Riemann sum approximated by an integral). This jump involves a 
logarithmic function with a cut along the linear segment [—iA, iA] . 

Now, to rigorously show that the approximation is valid one needs a very careful 
analysis of what happens near 0, because the "Riemann sum" approximation breaks 
down there. Indeed this is done in section 4.4.3 of [3], where, however, it is assumed 
that the loop (the jump contour) is located exactly at the S-curve (which solves 
the variational problem), and that this S-curve emanates from 0+ at an acute, 
non-right, angle to the horizontal axis. 

Even though in [3] the S-curve is assumed not to touch the linear segment 
[— iA, iA], this condition is actually only necessary locally near 0. There is no 
problem in allowing the S-curve to intersect the spike, as long as it emanates from 
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0+ at an acute, nonright, angle to the horizontal axis! Now it is proved in [3] that 
the S-curve solving the variational problem has to emanate from 0+ at an acute, 
non-right, angle to the horizontal axis. (This in fact follows immediately from the 
measure reality condition (4.18) applied at the origin.) So the proof of Theorem 
A.l will go through even if the S-curve intersects the linear segment [— iA, iA]. 
We then have the following. 

THEOREM 9. Consider the semiclassical limit (h — > 0) of the solution of (9)- 
(10) with bell-shaped initial data. Replace the initial data by the so-called soliton 
ensembles data (as introduced in [3]) defined by replacing the scattering data for 
tp(x,0) = tpo(x) by their WKB-approximation. Assume, for simplicity, that the 
spectral density of eigenvalues satisfies conditions (1). 

Then, asymptotically as h — > 0, the solution ip(x,t) admits a "finite genus de- 
scription", in the sense of Theorem A.l. 
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